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Abstract 



We present the result of a theoretical study of the quasi free electrofission of 
238 U. The exclusive differential cross sections for the quasi free scattering re- 
action stage have been calculated in PWIA, using a Macroscopic-Microscopic 
approach for the description of the proton bound states. The nuclear shape 
was parametrized in terms of cassinian ovaloids. The equilibrium deformation 
parameters have been calculated by minimizing the total nuclear energy. In 
the calculation the axially deformed Woods-Saxon single particle potential 
was used. The obtained single particle momentum distributions were aver- 
aged over the nuclear symmetry axis direction. The occupation numbers were 
calculated in the BCS approach. The fissility for the single hole excited states 
of the residual nucleus 237 Pa was calculated on the statistical theory grounds 
both without taking into account the pre-equilibrium emission of the particle, 
and with preequilibrium emission in the framework of the exciton model. 



I. INTRODUCTION 



Quasi free scattering of high energy electrons (QF) on nuclei is the field of nuclear physics 
which is traditionally devoted to the study of the single particle aspects of nuclear structure: 
single particle binding energies, momentum distributions, occupation numbers, etc [|l|. 

A new branch of these investigations is the study of decay channels of single hole states 
in the residual nucleus, created as a result of the QF process. Especially interesting is to 
study a fission decay following a QF process. In this case we have a single particle process in 
the first reaction stage, and essentially a collective process in the final reaction stage. The 
collective degrees of freedom are excited in the intermediate reaction stage due the residual 
interaction. 

This is a new sort of nuclear reaction which may allow one to get unique information on 
the dissociation of well defined single hole configurations ( which we can select by coincidence 
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(e, e'pj) into complex nuclear configurations, and its role in nuclear fission. The new and 
most important aspect of this reaction is that, after knocking-out a proton, we obtain the 
heavy nucleus 237 Pa in a single hole doorway state (see discussion below) which could undergo 
nuclear fission. Indeed, instead of dealing with collective doorway states, which are coherent 
sums of a great number of lp-lh configurations ( as the well-known giant resonances), these 
non-collective doorway states will be represented by only one, well defined, lh configuration. 
The residual interaction in 237 Pa mixes these lh configurations into more complicate 2h-lp 
and 3h-2p ones. So, there would be some competing channels for fission. It may occur 
either directly from lh configurations, or, with some delay, from mixed states (or their 
components). In a QF process we have in the initial state only one configuration; thus, the 
fission probability P/ should be more sensitive to the individual structure of this initial state 
as compared with conventional reactions, where the effects of the structure are averaged out 
over many single particle states forming the doorway. 

The unambiguous extraction of single hole contributions is possible only in an exclu- 
sive experimental scheme (reaction (e, e' pf)) and involves extremely thin targets (fission 
fragments have to leave the target with small energy losses), high energy resolution, and co- 
incidence requirement between the final particles in order to separate the single hole states. 
The exclusive (e, e'p/)-experiment is very difficult for practical realization, and never has 
been so far performed. The integral contribution of the quasi free electron scattering to the 
fission process was studied only in inclusive experiments: (e, /) p| and (e, e'f)) ||. These 
works dealt only with the issue of the QF contribution in electrofission. 

The advent of high energy, CW, electron accelerators combined with the development 
of high resolution facilities, opens the possibility of studying the fission channel for quasi 
free electron scattering in an exclusive experimental setup. The most suitable accelerator 
for this experiment is at the Thomas Jefferson National Accelerator Facility (TJNAF). 

For excitation of the residual nucleus to an well defined single hole state, the initial and 
final state interactions have to be negligible. This situation corresponds to high momentum 
transfer and high proton exit energy, when the Plane Wave Impulse Approximation (PWIA) 
for the calculation of the quasi free electron scattering cross section is valid. 

This work presents the results of a PWIA calculations for the quasi free (e, e'p)- 
differential cross section for deformed orbitals of 238 U, in the framework of the macroscopic- 
microscopic approach, plus an estimate of the fissility for single hole states in the residual 
nucleus 237 Pa, performed on the statistical theory grounds. These calculations could serve 
as first order magnitude guide-line for expected cross sections. 



II. PWIA CROSS SECTION 



In the first order Born approximation the electron with initial four-momentum = 
(ki,i£i) and final k 2fl = (A; 2, ^2) > transfers a virtual photon with four-momentum 
g M = (if, iuj) = fcx/i — k 2fJ , , resulting in the final state a knocked-out nucleon with 
Pn = (p£,iE p ) and a residual nucleus with Pa-u* — (Pa-Ij^-i)- 
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In the impulse approximation a virtual photon interacts with a bound nucleon (proton or 
neutron) of four- momentum p m = (~J? m , iE m ), which exits the nucleus with four-momentum 
Pfj. = (Ppj iEp) without further interaction (no FSI). The corresponding momentum diagram 
in the impulse approximation is shown in Fig. 1 for the Laboratory system.. 

In the plane wave impulse approximation (PWIA) jf m = — ~P a-i, and the missing 
quantities (momentum and energy of the proton before interaction) can be defined from the 
energy and momentum conservation law in the following way: 

~f m = l?p- (1) 
E m = io — T p — Ta-i, 

where E m = Ma~ \ + rn p — Ma is the proton missing (or separation) energy, T p is the kinetic 
energy of the outgoing proton, and Ta~i is the kinetic energy of the residual nucleus. The 
momentum and energy transfer of the virtual photon can be varied independently. 

In the PWIA the six folded differential cross section of the (e,e'p)— reaction in the 
Laboratory system has the following form 

d 6 a 

p p E p a ep S(E m ,p^), (2) 



dfl e dVt p ds2dE p 
where 

Vep = Vmott {V C W C + V T W T + VjWj + V S W S ) (3) 

is the off-shell electron-nucleon cross section, S{E m) p~^ l ) is the spectral function which 
defines the combined probability to find a bound proton with momentum p~^ on the shell 
with separation energy E m . 

The kinematic functions V in Eq. (^) can be expressed , neglecting the mass of the 
electron, as: 

V c = ^ (4) 
^ = ^ + tan 2 (|), (5) 

^ = |cos0^ + tan 2 (|), (6) 

^ = 4cos 2 + tan 2 (|), (7) 

and 

a 2 cos 2 % . 2sx . 2 6 e , , 
crmott = . 2 , 4 J l + — sin 2 -1 - 1 8 
4ef sin 4 ^ m p 2 

is the Mott cross section, is the angle between the scattering plane and the plane defined 
by the vectors "p* and if. 
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For the structure functions W in Eq. (|3]) we use the off-shell prescription of de Forest 



Wc = ^-{{E + E p f{F 2 + J^f) - g 2 (^i + ^ 2 ) 2 }, (9) 

W T = -2— (F 1 + k p F 2 ) 2 , 
2EE V ' " 



W z 



p 



+ E^F, + —K p F 2 
p^sin 2 7 / ql o 

where k p =1.793 is the anomalous magnetic moment of the proton in units of the Bohr 
magneton , 



E = \ PL + ml (10) 



m p is the mass of the proton, q = ("if, iu) , lj = E p — E, 7 is the angle between Jf p and 

Fi and F2 are the on-shell Dirac and Pauli proton form factors, respectively: 

2 

FiidD = —jr[GE(<ft + h^ G M{<l% (11) 

l + 4^f P 

K P E 2 (q) = —L-[G M (<Q - G E (q% (12) 



4m2 



where 

= (is) 

' 0.71 

^(^ = ^(9;)], (14) 

fi p = 2.793 is the proton magnetic moment in units of the Bohr magneton and q 2 in Eq.(|l~3"D 
is in (GeV/c) 2 . 

In the independent particle shell model the spectral function for the spherical orbitals 
a = nlj with binding energy E a takes the simple form: 

S(E m ,p^) =8{E- E a ) v 2 a n a (p^) , (15) 

where v 2 a and n a (p) are the occupation number and momentum distribution of the a orbital, 
respectively. The six folded (e, e'p)— cross section could be transformed into a five folded 
one: 

d 5 a 



dVL e dVt p dE p 



PpE p a ep v Q n a (pZ) , (16) 



where it is imposed energy and momentum conservation for the kinematics variables in a, 



Cff 
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III. SINGLE PARTICLE BOUND STATES 



The single particle bound state energies and momentum distributions were calculated in 
the framework of the macroscopic- microscopic approach by using the BARRIER code ||. 
The energy of the nucleus is presented as: 

E t ot = E LD + 5 E she u } (17) 

where E LD is the macroscopic liquid drop part of the energy and 5E she u is the shell cor- 
rection, which describes shell and pairing effects. Both shell correction and the macroscopic 
part of the energy have been calculated according to |J. 



A. Nuclear shape parametrization 



Only axially symmetric nuclear shapes have been considered in the present work, and 
the deformed shape (up to and beyond its separation into two fragments) can be conve- 
niently described by the Cassini ovoids [BUil. The potential-energy surfaces are calculated as 
functions of e (elongation) and (hexadecapolar momentum). From these potential energy 
surfaces, the equilibrium (ground state) deformation parameters e and were calculated 
by minimizing the total nuclear energy ( Eq.flTTD): e = 0.227 and a 4 = 0.059. 



B. Nuclear Potential 

An Woods- Saxon potential f|, consisting of the central part V, spin-orbit Vso, and the 
Coulomb potential Vc ou i for protons, was employed: 

V ws (r, z, e, a) = V(r, z, e, a) + V so (r, z, e, a) + V Cou i(r, z, e, a) (18) 

The real potential V(r, z,e,a) involves the parameters V , r and a , describing the depth 
of the central potential, the radius and the diffuseness parameter, respectively, and it is 
expressed as: 

V{r, z, e, a) = — =-, (19) 



1 + exp 



dist(r,z,e,a) 



where dist(r, z,e,a) is the distance between a point and the nuclear surface, and e and a 
are deformation parameters. 

The depth of the central potential is parametrized as 

V = V [l ± k(N - Z)/(N + Z)], (20) 

with the plus sign for protons and the minus sign for neutrons, with the constant k = 0.63. 
The spin-orbit interaction is then given by: 

V so (r,z,e,a) = A (^^j W(r,z,e,a) ■ (a x p), (21) 
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where A denotes the strength of the spin-orbit potential and M is the nucleon mass. The 
vector operator a stands for Pauli matrices and p is the linear momentum operator. 

The Coulomb potential is assumed to be that corresponding to the nuclear charge (Z — 
l)e, taken to be uniformly distributed inside the nucleus. It is computed in cylindrical 
coordinates by using the expression given in 0. 

C. Single particle potential parameter definitions 

For the ground state deformation of 238 U, small changes in A (spin-orbit potential 
strength) and r _ so (spin-orbit potential radius) of the Chepurnov parameters [f| are in- 
troduced in order to reproduce adequately the spin/parity of the levels sequence. Using 
single particle states obtained by this procedure, the quasiparticle states can be calculated 
for the first minimum region, providing spin, parity, energy and level spacing for the ground 
and some low-lying states. The quasiparticle spectrum was obtained by using the semi- 
microscopic combined method JT0[ . 



The potential parameters were chosen to give the best fit to the spectrum of single- 
quasiparticle excitations of the Z-odd neighboring nuclei 239 Np. 



D. Single particle wave functions 

The Hamiltonian matrix elements are calculated with the wave functions of a deformed 
axially symmetric oscillator potential. The wave functions in the coordinate space 0, are 
expanded into eigenfunctions of the axially deformed harmonic oscillator potential. 

These eigenfunctions form a complete orthonormal basis for the single particle Woods- 
Saxon wave function 

*,(3,a)= £ C; p ^ AE $ n ^ AE (it a) . (22) 

n p ,n z ,A,E 

From this expansion, we may conveniently express the single particle Woods-Saxon wave 
function in momentum space: 

%(t,a)= C n p ,n z A,^n p ,n^(t,a), (23) 



n p ,n z ,A,E 



with 



*,v,..al ( II. a) = / dle-'K X^n^ ( R.cr) (21) 



normalized to one. 



We define densities n« (T? j in momentum space in an analogous way of that in the 
configuration space: 







\®t(r,z) 2 + 


®i (r,z)\ 



(25) 
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with 



and 



with 



$f (r, z) 



n„,n z ,k,T, 



n 



(if) = m (k, k z ) = |$+ (A;, A;,) | 2 + |$r (fc, fc 2 ) 



(26) 



(27) 



*?(*,*.) = 4= ^^W^^aA^a^P^)- (28) 

np,n z ,A,2j 

These single particle momentum distributions rij (j^^j were averaged over nuclear sym- 
metry axis directions. 

Similarly to the total density 



p = E 2w ? * (^) > 

i 

the total momentum distribution is given by 



(29) 



(30) 



where t> 2 is the occupation probability resulting from the BCS model. [n|, JTD 
The results for the occupation number calculations are shown in Fig. 2. 
The energies of the 238 U proton bound states are shown in table 1. 



IV. FISSILITY 



The quasifree knockout of nucleons leads to the excitation of the residual nucleus. This 
excitation energy (E*, nucleus A-l) has two origins: holes in the shells of the nucleus A, 
which appear as a result of the knockout of nucleons, and final state interaction (FSI) of the 
outgoing nucleon, which we assume as negligeable due to the high energy of the proton. 

The fast, quasi free reaction stage occurs at zero thermal excitation (ground state) of 
the initial nucleus 238 U, and results in a single hole in one of the shells. This single hole 
configuration forms a doorway for a thermalization process which leads to the thermal 
excitation E*of the residual nucleus 237 Pa. 

The thermalization is a complicate process which involves creation of new many particle- 
hole configurations in competition with particles emission and fission, and for some doorway 
configurations it might has non statistical character, but, as a first guide-line for order 
of magnitude estimates we calculate the total fission probability (nucleus with energy E* 
deexcites in several steps) on the statistical theory grounds, both with and without taking 
into account the preequilibrium decay. 
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A. Compound nucleus model 



Firstly, we considered a extreme situation, by assuming that the residual interaction 
leads to thermalization and formation of compound nucleus just after the fast reaction 
stage, without any preequilibrium particle emission In this case, the compound nucleus 
excitation energy is assumed to be : 



E* 



-E n 



(31) 



where E a is the energy of the bound state (hole). 

For calculations of compound nucleus fissility we used the Bohr Wheeler [O] and 



Weisskopf |T3[ models for the description of the evaporation/fission competition. It was 
developed a Monte Carlo algorithm for the evaporation/fission processes which includes not 
only the neutron evaporation vs. fission competition, but also takes into account the proton 
and alpha-particle contributions. 

The probability for the emission of a particle j with kinetic energy between E k and 
E k + dE k is calculated within the Weisskopf statistical model [T^| as: 



P f (E k )dE k = 1] a j E k dE k 



(32) 



where <7j is the nuclear capture cross section for the particle j, = -Tp, where g denotes 
the number of spin states, and m is the particle mass. The level densities for the initial and 
final nucleus, pi and pf, respectively, are calculated from the Fermi gas expression 

p(E*) = exp [2( aE*) 1 ' 2 ] , 

where a is the level density parameter (see below), 



E* 



E* - (B< + Vn 



(33) 



E* is the nuclear excitation energy in the initial state, Bj is the particle separation energy, 
and Vj is the Coulomb barrier corrected for the nuclear temperature, r, defined by E* = ar 2 . 
The particle emission width is calculated as 



Pj{E k )dE k . 



(34) 



From this general equation, the fc-particle emission probability relative to the j-particle 
emission is: 



7fc E* k a i 

lj E j a k 



exp 



2 Ua k Elf 2 - ( aj E*) 



1/2 



(35) 



The level density parameter for neutron emission is [ H |: 

a n = (0.134 A - 1.21 ■ 10~ 4 A 2 ) MeV" 1 , 



(36) 



S 



and for all other particle emission this quantity is related to a n by, 



Clj TjCl n) 



(37) 



where Tj is an adimensional constant. 

Shell model corrections are not taken into account. For high excitation energies their 
effects are likely to cancel each other upon averaging over all possible nuclei created during 
the reaction. 



from Weisskopt |0, we get 



Using the fission width from the liquid drop model [ 12| , and the neutron emission width 

(38) 



K, exp [2 ((a, E}) 1/2 - (a n El) 1 ' 2 ) 



where 



K f = K a n - 



2((a f E})^ - 1 



4A 2 / 3 a f E. 



(39) 



with K = 14.39 MeV and E* = E* — Bf. Here Bf is the fission barrier height discussed 
below. 

For proton emission we get 



r„ 



e; 

E* 



exp 



2(a„) 1/2 (M;) 1/2 -(K) 1 / 2 ) 



(40) 



and for alpha-particle emission |TB| 



r, 



2 E* 



E* 



f exp [2(a n ) 1 / 2 ((r a K) 1 / 2 - (K) 1/2 )] 



(41) 



In the above equations, the Coulomb potential for protons is [[17 



and for alphas, 



v - c k A^3l m 

Vp ~ r (A-4)V3 + jR/ ^ 



Va ~ r (A- 4)V3 + i? a ' (43) 



where K p = 0.70 and K a = 0.83 are the Coulomb barrier penetrability for protons and alpha 
particles, respectively, R p = 1.14 fm is the proton radius ,R a = 2.16 fm is the alpha particle 
radius, and ro = 1.2 fm . The factor C introduces in a semi-empirical way the dynamical 



effects in particle separation energy and fission barrier due to the nuclear temperature Hl7 
namely 

VP* 

C = 1 -^> ( 44 ) 
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where B is the total nuclear binding energy (B=1794 MeV for 237 Pa ]|i~7| ' 



The fission barrier is calculated by [17 



B f = C(0.22(A - Z) - 1.40Z + 101.5) MeV. 



(45) 



The neutron separation energy was taken as 5.78 MeV for the first step ( Pa) , and 



for the other steps as [fLq| : 

B n = (-0.16(4 -Z) + 0.25Z + 5.6) MeV, (46) 
while the proton and alpha-particle separation energies are calculated through the nuclear 



mass formula 18 



B p = m p + M(A- 1,Z- 1) -M(A,Z), (47) 

where m p is the proton mass, and M(A, Z) is the nuclear mass calculated with the param- 
eters from reference [18]. For the alpha particles we get 

B a = m a + M(A-A,Z -2)- M(A,Z), (48) 

where m a is the alpha particle mass. 

These values reproduce the experimental data for Pf (see discussion below) 
The present Monte Carlo code for evaporation-fission calculates, at each step % of the 
evaporation chain, the fission probability, F i: defined as 

77 1 v 1 n 



' i + (Et), + (Jk) i + (J*),' 

An evaporating particle j is randomly chosen (neutron, proton or alpha particle), 
according to its relative branching ratios. Once one of these particles is chosen , the mass 
and atomic numbers are recalculated through 

A i+1 — Ai — A A,, 

and 

Zi+i = Zi — AZi, 

where AAi, and AZj, are, respectively, the mass and atomic numbers of the ejected 
particle at the ith step in the evaporation process. The nuclear excitation energy is modified 
according to the expression 

where Bi and Tj are the separation and the asymptotic kinetic energies of the particle being 
ejected, respectively. For neutrons T = 2 MeV, and for protons and alpha particles T = 
MeV. The expressions described above ensure that the nuclear excitation energy will be, 
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at each step in the evaporation chain, smaller than in the previous step. This process 
continues until the excitation energy available in the nucleus is not enough to emit any one 
of the possible evaporating particles. At this point the evaporation process stops, and we 
can calculate the nuclear fissility by the expression 



.i=o 



(49) 



Using the model described above, we calculated the fissility for 237 Pa (figure 3, solid 
curve). Peaks observed for the fissility reflect the opening of the fission channel in the 
daughter nuclei. Figure 3 also shows ( rectangle) the data for the fissility of 237 Pa obtained 
by extrapolation of the neutron to fission widths ratios for Z = 91, and A = 230, 231, 232, 233 
to A = 237 ||19|| , by using the empirical trend presented in Vandenbosch and Huizenga (16 



It should be pointed out that in our calculations of the fissility we assumed that the 
hole excitation energies for an A-l nucleus correspond to the compound nucleus excitation 
energies , that is to say, the complete thermalization is reached without any preequilibrium 
decay. Such calculations could be considered as an upper limit estimate for the fissility. 



B. Exciton model 



During the thermalization of the hole excitation energy, the nucleus A — 1 could undergo 



particle evaporation (preequilibrium decay [pG|| , plj ). 



In this case, the energy of the hole is not attributed to the nuclear temperature but, 
instead, assumed as a characteristic of the doorway state in the thermalization process 
followed by the emission of particles or fission. 

The calculation involving the preequilibrium decay was performed within the framework 
of the exciton model using the code STAPRE. In this model, the states of the system 
are classified according to the number of excitons n, which corresponds to the total number 
of excited particle p and hole h degrees of freedom, n = p + h. The exciton model included 
in STAPRE does not distinguish between protons and neutrons. Starting from a simple 
configuration of low exciton number, the system is assumed to equilibrate through a series 
of two-body collisions and to emit particles from all intermediate states. The application of a 
two-body interaction to the states of a (p, h) configuration results in states with (p+1, h+1), 
(p,h), and (p — l,h — 1) excited particles and holes. The difference between the number 
of excited particles and holes remains fixed, justifying the use of the exciton number to 
label the states. However, the transition rates, which are an averaging over all states of 
a configuration, do depend on the number of excited particles and holes. The equation 
governing the time development of the occupation P(n) of the n — th exciton configuration 
can thus be written as 

= A_(n + 2) P(n + 2) + A (n) P(n) + X + (n - 2) P(n - 2) - X(n) P(n) (50) 

where A(n) is the total transition rate, 
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A(n) = A_(n) + A (n) + A+(n) + A e (n), 



(51) 



with A e (n) being the total rate of particle emission from the n — th exciton configuration. 
The quantities A_(n), Xq(ji), and A + (n) are the average rates for internal transitions from 
the n — th exciton configuration with a change of exciton numbers by -2, 0, or +2. 

The internal transition rates can be written as the product of the average squared matrix 
element of the residual interaction |M| 2 with the relative density of available states. For the 
latter, STAPRE uses the expressions of Williams |E3| as corrected for the Pauli principle by 



Cline [p4 |. These yield 

A» es \_(p,h,E) = ^ W 9(9E-C^ +1 )^ 

n p + h + 1 

X (n) = X (p,h,E) = ^\M\ 2 g(p + h-l)(gE- C p , h ) , 
A + (n) = A+(p, /i, £?) = — |M| 2 OT % + h - 2), 

where 

C P ,u = \{P 2 + h 2 ) , 

with E the excitation energy of the system. The parameter g is the single-particle state 
density, which is taken to be g = -%a, with a as the level density parameter. Following 



Kalbach-Cline [2S], the average matrix element is approximated as 



\M[ 



A 3 E : 



where A is the mass number of the system and fu is a parameter, which we assumed to be 
fu = 230 MeV 3 in our calculations. 

The particle emission rate X e (n) is the sum of the integrated proton and neutron differ- 
ential emission rates, X eu (n, e) de u , which are determined through considerations of detailed 
balance |2lj , 



X ev (n, e) de u = \ eu (p, h, E, e v ) de v (53) 
1 t \ d u{p-l,h,E - B v - s„) 

\\^v &v ®v\£v) ti-v 7 7 — ae vi 



n 2 h^ v " uy v ' v u(p,h,E) 

where \i v is the reduced mass of the emitted neutron/proton, e v its outgoing kinetic energy, 
B v its separation energy, and a v {e u ) is the cross section for the inverse absorption process. 
The factor R u is a simple correction standing for the fact that neutrons and protons have 
not been distinguished in the process; thus, 



J N/A for neutron emission 
1 Z/A for proton emission 



The densities of states are taken to be the Williams densities, 
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(54) 



where the Pauli blocking correction is 

A V)h = i (p 2 + h 2 + p - 3h) . 

The differential emission rates differ from those of usual Weisskopf compound nucleus emis- 
sion by the factor R v and by the use of exciton state densities rather than compound nucleus 
ones. 

The time evolution equation, Eq. (|5CiD , form a set of coupled linear differential equations 
whose solution could be written in the form of a vector as, 

Pit) = exp[-At] P , 

where the matrix A is given by 

A nn / = \{n) 8 n > >n - A_(n + 2) 5 n ^ n+2 - A (n) 5 n >, n - X + (n - 2) 5 n / jn _ 2 , 

and the vector P describes the initial exciton configuration of the system, 

P (n) = P Q (p,h) = 5 PtP0 5 hM . 

The decay of the system into all possible final configurations can be obtained by integrating 
the total emission rate over all time, 

P CO 

J2 \ e (n)P(n,t)dt = J2Un) (A" 1 )^, P (n'). (55) 

n ^° n,n' 

The decay of the fraction of the initial probability, which survives preequilibrium emission, is 
described using the Hauser-Feshbach formalism. We have considered fission in competition 
with neutron and gamma emission. 

The initial configuration in 237 Pa consists of one-particle at the Fermi level and one-hole 
in a bound state. This configuration is consistent with the proton knockout reaction for 
238jj initiating the statistical cascade. Our calculations were performed assuming an one 
-hole initial configuration of the 1=0 partial wave alone. The particle at the Fermi level 
contributes negligibly to the equilibration process. The fission barriers, neutron separation 
energies and level density parameters were taken the same as those of the compound 
nucleus calculations in the previous section. 

The exciton model fissility results for single hole states of 237 Pa are shown in fig 3 by the 
dotted curve. We note that these calculations for fissility show a smoother behavior than 
that for compound model. The preequilibrium particle emission removes some excitation 
energy before a equilibrium is reached reducing, therefore, the probability of opening new 
chances for fission. 
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V. FINAL RESULTS 

The differential cross section for the (e, e'p/) -react ion was obtained by assuming an 
isotropic angular distribution for the fission fragments, and the fissility as a factor : 

*° 1 *" -P, (56) 



dfl P dVtr, dEr, dQ f 4n dVL e dQ„ dE„ 



Fig. 4 shows the seven folded differential cross sections dn dn dE dn ^ for some bound 
proton states (table 1) and the compound nucleus model fissility (solid curve in fig. 3) calcu- 
lated for e 1 =2000 MeV , 8e = 23° and the parallel kinematics j27[. In this kinematics ex 
and 6e are fixed, and for each value of e 2 the proton spectrum is measured in the direction 
of if , varying each time the angle 6^^. (see fig. 1). For such scheme of measurements the 
initial (missing ) momentum of the proton is always parallel (or antiparallel) to if . 
The parallel kinematics simplify the accounting of FSI, since there are no contribution of 
interference terms in the cross sections (see Eqs. ([|)). Figs 5 and 6 show the momentum 
distributions for states used in the calculation of the cross sections presented in fig. 4, and 
fig. 7 shows the outgoing proton kinetic energy and angle Qp p k Y versus e 2 for the parallel 
kinematics we use. 

Figs. 8 and 9 show the differential (e, e'pf)— cross sections calculated for the same E\ 
and 9e but for two fixed proton angles: Qp v k x = 0.98 rad for the group of proton states of 
238 JJ which have a maximum in the low missing momentum region (fig.5), and Qp v k x = 0.82 
rad for the group having a maximum in the high missing momentum region (fig. 6) . These 
angles dp^were chosen in order to achieve parallel kinematics, that is 6p q ~ , and 
maximum for cross sections at both the low ( 9p = 0.98 rad ) and high (dp = 0.82 rad ) 
missing momentum regions. Fig. 10 shows the missing momentum P m and the angle 9p pQ as 
functions of the outgoing proton kinetic energy for 6 P = 0.98 rad ( solid curve) and 0.83 rad 
( dashed curve). 

It is seen from the figures 8,9 and 10 that for such a choice the cross sections have 
maxima at E p around 300 and 400 MeV and, for these energies, the proton angles 8p pq are 
small (parallel kinematics). 

The differential cross-sections presented in figures 4, 8 and 9 correspond to the situation 
when the hole excitation energies for an A-l nucleus is the compound nucleus excitation 
energies , that is to say, the complete thermalization is reached without any preequilibrium 
decay. Such calculations could be considered as an upper limit estimate for the cross section 



VI. CONCLUSIONS 

We presented a theoretical study for the quasifree electrofission of 238 U. 

The proton bound states were calculated in the framework of the Macroscopic- 
Microscopic approach, using the axially deformed Woods-Saxon single particle potential. 
The occupation numbers were calculated in the BCS approach. 
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The exclusive differential cross sections for the quasi free scattering reaction stage were 
calculated in PWIA using off-shell electron-nucleon cross sections. 

The fissility for the single hole states of the residual nucleus 237 Pa was calculated in the 
framework of two approaches: compound nucleus model without taking into account the 
preequilibrium emission of the particles, and the exciton model accounting for preequilibrium 
emission. Both models exhibit the same general trend, but the fissility as given by the 
preequilibrium model is smoother. 

The obtained results could serve as a first guide-line on order of magnitude estimates of 
the expected cross sections for quasi free electrofission of 238 U. 
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Table I: Proton single-particle levels of 238 U. The Fermi level is the level 46. 
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Figure Captions 

Fig.l. Momentum diagram of the A(e,e'p)A — 1 reaction in the impulse approximation: 
k i, and k 2 are the initial and final electron momenta, respectively; ~p* m is the momentum 
of proton before interaction; is the momentum of knocked-out proton. 
Fig. 2. Occupation probabilities for the single particle bound states of 238 U. 
Fig. 3. Fissility of 237 Pa vs the hole excitation energy. The solid curve shows the compound 
nucleus model calculation, assuming that the hole excitation energies correspond to the 
compound nucleus excitation energies: the complete thermalization is reached without any 
preequilibrium decay. The dotted curve corresponds to the exciton model calculations, 
which take into account the preequilibrium decay. The rectangle shows the extrapolated 
experimental data (se text for details). 

Fig. 4. Seven folded differential cross sections for the parallel kinematics. The calculations 
of the cross sections were accomplished for E\ = 2000 MeV, 6e = 23°, = 0°. The cross 
section for E=-10.388 MeV state ( subbarier) is multiplied by 100. 

Fig. 5. Momentum distributions for some proton bound states having maxima at the low 
missing momentum region. 

Fig. 6. Momentum distributions for some proton bound states having maxima at the high 
missing momentum region. 

Fig. 7. Variation of the angle Ppkl (see fig.l) and the outgoing proton kinetic energy E p 
vs e 2 for the parallel kinematics and e\ = 2000 MeV, 6e = 23°. 

Fig. 8. Differential cross sections for some bound states of 238 U, having maxima of the 
momentum distributions at the low missing momentum region. The calculations of the cross 
sections were accomplished for E\ = 2000 MeV, 9e = 23°, = 0° and 9 Ppkl = 0.98 rad. 
Fig. 9. The same as in figure. 5, but for some bound states of 238 U, having maxima of the 
momentum distributions at the high missing momentum region and for Ppkl = 0.82 rad . 
Fig. 10. Missing momentum p m and angle 9 Ppq as functions of the outgoing proton kinetic 
energy for 9 Ppkl = 0.98 rad ( solid curve) and 0.83 rad ( dashed curve). 
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